Structural crossover of polymers in disordered media 
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We present a unified scaling theory for the structural behavior of polymers embedded in a disor- 
dered energy substrate. An optimal polymer configuration is defined as the polymer configuration 
that minimizes the sum of interacting energies between the monomers and the substrate. The fractal 
dimension of the optimal polymer in the limit of strong disorder (SD) was found earher to be larger 
than the fractal dimension in weak disorder (WD). We introduce a scaling theory for the crossover 
between the WD and SD limits. For polymers of various sizes in the same disordered substrate 
we show that polymers with a small number of monomers, N <^ N* , will behave as in SD, while 
large polymers with length N ^ N* will behave as in WD. This implies that small polymers will 
be relatively more compact compared to large polymers even in the same substrate. The crossover 
length A^* is a function of u and a, where v is the percolation correlation length exponent and a is 
the parameter which controls the broadness of the disorder. Furthermore, our results show that the 
crossover between the strong and weak disorder limits can be seen even within the same polymer 
configuration. If one focuses on a segment of size n <g N* within a long polymer {N ^ N*) that 
segment will have a higher fractal dimension compared to a segment of size n 3> A*. 



The study of polymers [H, S, H, i, i, i, 0, El in the 
presence of disordered media is of broad scientific interest 
and is relevant to many fields such as protein folding [ij, 
polymers in porous media ^ and spin glasses 0] . It also 
has many important applications such as in enhanced oil 
recovery, drug delivery and DNA sorting A polymer 
embedded in a disordered energy substrate, at low tem- 
peratures, will settle down in the optimal configuration, 
i.e. the configuration with the minimum energy. 

In this paper we study the scaling of the structural 
properties of such an optimal configuration. A linear 
polymer of N monomers, under disorder, can be modeled 
by an N step self avoiding walk (SAW) on a lattice where 
each site (or bond) in the lattice is assigned an energy e, 
taken from a given distribution. The optimal configura- 
tion in such a model is the SAW of length N for which 
the sum of the energies along its path is minimal. For 
generating a broad disorder it is common to use the dis- 
tribution P(e) — l/(ae) (e < e°) where the parameter a 
controls the broadness of the disorder 12, 13]. Such 
a distribution generates an exponential disorder where 
the energy value of each site i on the lattice is given by 
e-i = exp(ari) where is a random number taken from 
a uniform distribution between [0, 1]. For a — > oo we ob- 
tain the strong disorder (SD) limit and for small values 
of a the weak disorder (WD) limit. In the WD limit es- 
sentially all sites contribute to the total sum of e^, while 
in the SD limit the total sum is dominated by a single 
site with the maximum energy [^, [lo| . 
Smaller et. al. [ll| studied the properties of a linear 
polymer in the WD limit using both a uniform and a 
Gaussian energy distribution. They found that for 2D 
the end-to-end distance of the polymer, R, scales with N 
as iV ~ R'^^opt, where the fractal dimension is d^^j ~ 1.25. 



This result is different from that found by Braunstein 
et. al. [l^ for the SD limit. Braunstein et. al. used 
the exponential disorder (e^ = exp(ari)) and obtained 
N ~ Ji<iopt -^vith dopt — 1-5. Thus, the polymer is more 
compact in the presence of strong disorder compared to 
weak disorder. 

In this paper we present a unified scaling theory for 
the crossover between WD and SD limit. We claim that 
the crossover between strong and weak disorder for the 
optimal polymer problem depends on the characteristic 
size a'^ where v is the correlation length exponent from 
percolation theory 1^, 3, 13] ■ 

This claim can be explained as follows: Consider an 
infinite lattice where each site i has a value pi where pi G 
Jo, 1]. In a percolation process an increasing fraction p of 
the sites with the lowest values, are occupied (the others 
are removed) until the point p = Pc is reached where 
the lattice undergoes a phase transition. For p > pc an 
infinite number of the lattice sites are connected while for 
p < Pc the lattice is separated into small finite clusters. 

Now assume each site is associated with an energy 
ei — exp{api) and an infinite polymer is observed on the 
lattice. In SD the highest energy, tmax , along a polymer 
configuration dominates the total energy. The energy 
£max which is the highest energy along the optimal poly- 
mer configuration, cannot be smaller than exp{apc) and 
cannot be larger than exp{apc) and therefore must be 
equal to exp{apc). If imax < exp{apc) then according to 
percolation theory only a cluster with a finite number of 
sites exists and an infinite polymer cannot be embedded 
inside this cluster. If emax > exp{apc) then according 
to percolation theory the polymer is not optimal since 
there exists a polymer with a lower energy configuration 
for which emax = exp{apc). Thus, the optimal configura- 
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FIG. 1: Simulation results for the scaling of /?iv(a)/7?Ar(oo) 
as function of N/a'^ for different values of o. The dashed 
lines represent the asymptotic expected value of the slope 



tion is achieved for (.max = exp{apc) since pc is the lowest 
value for which an infinite cluster that all its sites have 
energies lower than Cmax exists. 

For a finite lattice of size R, the value pc for which the 
percolation transition occurs is distributed as a narrow 
distribution with a mean pc and a standard deviation 
of (T = R^^^" Let ei — exp{api) be the largest 

energy and £2 = exp{ap2) the second largest energy on 
the lattice. Since log{^) ~ a{pi — P2) and pi and p2 are 
both from the same narrow distribution it follows that 
logi^) ^ a{pi -P2) ^ aa ^ aR'^/" . Therefore, the SD 
limit (ei/e2 S> 1) is obtained for R <^ a" while the WD 
limit is obtained for R';^ . 



Similar 
Buldyrev et. 



considerations 
al. 



have been presented by 



14l | for the crossover from SD to WD 
in the optimal path problem. Yet, the optimal polymer 
problem is significantly different from the optimal path 
problem. The optimal path problem explores the optimal 
path between two fixed points (fixed R) where the path 
length varies. In the optimal polymer problem, the poly- 
mer length is fixed and one tries to minimize the sum of 
the interaction energies between the monomers and the 
substrate, while the end-to-end distance varies. While 
the complexity of the optimal path problem is 0{N'^), 
finding the optimal polymer is an NP (Non-deterministic 
Polynomial time) problem since nearly all the possible 
configurations need to be explored in order to find the 
optimal configuration. To reduce the computational time 
of our simulations we have used the following optimiza- 
tions: a) During the calculation of the optimal polymer 
of length TV, the optimal polymers of length 1 to are 
also calculated, b) A new site in the lattice is explored 
only if the path until that point does not exceed the op- 
timal total energy of length N obtained until that point. 
The second optimization is only useful for the SD limit, 
where one energy dominates the total sum 18|. On the 



other hand for the SD limit the accuracy of the float- 
ing point (double) in conventional computers is limited 
to a 36. We have overcome this limitation by using 
external software that enables unlimited floating point 
accuracy, but demands signiflcantly more computational 
time. These computational problems are the reasons why 
our simulation results are limited to polymer conflgura- 
tions of length TV = 50. 
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FIG. 2: Simulation results of _Rjv(a)/_R]v(oo) as function of 
TV for several values of a. For small values of TV all the poly- 
mers behave as in the SD limit and independent of the dis- 
order have the same fractal dimension. For larger values of 
TV, i?iv(a)/_Rjv(oo) increases towards the WD region. The re- 
sults are limited to TV = 50 because calculating R{a) for large 
a slows down he calculations significantly 

Next, we propose a single scaling function for the de- 
pendence of i? = -Rat (a) on TV that includes both the 
SD and WD regions. We expect that the scaling func- 
tion will depend on the ratio R/a'^ which represents the 
relative strength of the disorder. Since the polymers 
are of length TV we use instead of R/a"^ the scaling pa- 
rameter TV/a" = N/a'"^""* which is related to R/a^ by 



Thus, we propose that 



TV 

-Rjv(a) ~ Rn {00) f ( —7 



with ,f{x) a scaling function given by 



const. X 



X 



(1) 



(2) 



where d. 



opt 



1.5 (SD) , d\ 



opt 



1.25 (WD) and 

i?Ar(00) - TVl/''""'. 

Indeed from Eq. H]) and ^ follows that for TV < a"' 
(SD) Rn N^/'^-p^ while for TV > a'" (WD) hmit we 
obtain Rn ~ 7V^/'^°p* the known result for the WD limit. 
Fig. [1] presents simulation results supporting the scaling 
function of Eq. ^ and ([2]). The SD region can be seen 
clearly in the flat area for N/a'^ <C 1, while the sharp 
increase in the slope indicates the crossover to WD until 
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the point where the slope fits the asymptotic expected 
value, l/d'gpt — ^/dopt (clashed hne), see Eq. ([2])- 

Until now we have been relating to SD and WD as if 
the polymers were in different substrates. But actually 
the crossover between the SD and WD is present even 
on the same energy substrate as can be seen from the 
scaling function of Eq. ([T]) and ([2]). In other words, 
for a given a, a small enough polymer (A'^ <^ ) will 
behave as in SD, while a large enough polymer {N :s> a" ) 
will behave as in WD. This property is demonstrated 
in Fig. [5] where we plot the ratio i?Ar(a)/i?Ar(oo) as a 
function of N for different values of a. For small values of 
N the polymers behave as in SD since RN{a) / Rn{oo) ~ 1 
(recall that Rpf{oo) is the SD limit). In contrast, for 
iV ^ 1 a crossover towards WD occurs represented by 
the larger R values. Note, that the crossover value N*{a) 
increases for larger a. 
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FIG. 4: Simulation results for the scaling of Rn{a) / Rn{oo) as 
function of n/a'' for different values of a. The different values 
of n represent different segments within the same polymer. 
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FIG. 3: (a) Illustration for the crossover between WD and SD 
inside the same polymer. Inside a segment of length n < a" 
the polymer is more compact (has a higher fractal dimension) 
compared to a larger segment or the full polymer. 



FIG. 5: Simulation results for the cumulative probability den- 
sity function, p{R, N, a) as function of R for different values 
of iV and a in the region 0.075 < N/a"' < 2. 



Next we show that such a crossover also exists within 
the same polymer configuration. We argue that a seg- 
ment of monomers of length n <C a'^ , within a polymer 
of length N ^ a'^ , will behave as in SD and will be more 
compact compared to a segment of size n' 3> a'^ that will 
behave as in the WD limit. An illustration of this prop- 
erty is shown in Fig. [31 In Fig. [3] we present simulation 
results showing that the same scaling relations describing 
the crossover between SD and WD for polymers of differ- 
ent sizes is also correct within the same polymer. Here 
the simulations were performed for different segments of 
size n of the same polymer and Rn{o.) is the end-to-end 
distance of these segments. 

When studying the probability density function of 
polymers p = p{R,N,a) we obtain an interesting result: 
The width of the distribution does not depend on the 
crossover parameter of the disorder N/a'^ , only on the 
normalization parameter R/{R). This is different from 
the results obtained in optimal path where the proba- 
bility density function depends on the disorder crossover 
parameter [14]. This difference can be understood by the 



fact that the polymer configuration is of fixed length N. 
While the length of an optimal path between any two 
nodes within a distance R may be very long, enabling a 
broader distribution of the path lengths in SD, the end- 
to-end distance R in the optimal polymer configuration is 
limited from above by the fixed length N (in optimal path 
one measures the distribution of the length, N, while in 
polymers we measure the distribution of R). Therefore, 
the scaling relation for the probability density function of 
the optimal polymer depends only on the normalization 
parameter R/{R) while the dependence of iV and a are 
determined by (i?): 

p{R,N,a) = fij^) (3) 

Indeed, Fig. [5] shows simulation results for the cumu- 
lative distribution of p{R, N, a) for different values of 
-R, a that correspond to both the strong and weak 
disorder regions. In Fig. [S] the scaling of these lines is 
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FIG. 6: Simulation results for the cumulative probability den- 
sity function, p{R/{R),N,a) as a function of R/{R). When 
normalizing by R/{R) the different curves collapse into a sin- 
gle curve. 



presented where all the different distributions both for 
SD and WD collapse into a single curve. 

Summarizing, we have presented a unified scaling the- 
ory for the optimal polymer configuration, embedded in a 
disordered energy substrate. The structure of the poly- 
mer at the minimal energy depends on the strength of 
the disorder. We find that the crossover between the SD 
limit and the WD limit occurs at N*{a) = N/a"' . For 
N/a"' > 1 the WD is obtained while for N/a"' < 1 
the SD limit is obtained. We present simulation results 
showing that this transition occurs even inside the same 
polymer. Therefore, a segment of n monomers where 
n/a'^ 1 (the SD limit region) will be more compact 
compared to the full polymer where N/ a'^ ^ 1 (the WD 
region). 

Our results of the crossover between the SD limit and 
the WD limit may be observed in neutron scattering ex- 
periments. Since the structure factor decays with a power 
(equal to the fractal dimension) of the wave vector, we 
expect to see a crossover in the structure factor at the 
wave length for which the transition between the SD and 
WD occurs. For 3D. the crossover would occur from 



^opt — ^ 



1.82 in SD limit hj to d' 



opt 



1.4 in the WD 



limit [11]. Recently, experimental studies of the scaling 
properties of the linear chain configuration of DNA knots 
adsorbed onto a mica surface under weak trapping (20| 
have shown that R scales with N with a fractal dimen- 



sion of 1.51. This value is significantly different from the 
classical Flory result df = 4/3. The fractal dimension 
df = 1.51 found in 2D might be explained by Eq. ([T|) in 
the limit of strong disorder. 
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